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We formulate statistical mechanics for a mixed diamond chain with spins 1 and 1/2. 
Owing to a series of conservation laws, any eigenstate of this system is decomposed into 
eigenstates of finite odd-length spin-1 chains. The ground state undergoes five quantum phase 
transitions with varying A, a parameter that controls frustration. Wc obtain the values of 
the residual entropy and Curie constant which characterize each phase and phase boundary 
at low temperatures. We further find various characteristic finite-temperature properties 
such as the nonmonotonic temperature dependence of magnetic susceptibility, the multipeak 
structure in the A-dcpendcncc of entropy, the plateau-like temperature dependence of entropy 
and the multipeak structure of specific heat. 

KEYWORDS: mixed spin chain, diamond chain, frustration, thermodynamics, residual entropy, 
exact solution 

1. Introduction 

The quantum fluctuation effects of one-dimensional quantum magnets have been exten- 
sively studied in various geometrically frustrated lattices. Owing to the interplay between 
quantum fluctuation and frustration, various exotic quantum phenomena can take place. 
For example, spontaneous dimerization,^) 1/3-plateau with spontaneous trimerization,^"^) the 
transition between quantum and classical plateaus, singlet cluster solid states,^-* quantized 
and partial ferrimagnetisms^"^^ and spin quadrupolar phases^' have been reported. 

Among a variety of models and materials with strong frustration, the uniform diamond 
chain (UDC) has been attracting the interest of many condensed matter physicists. The UDC 
consists of successive diamond-shaped units, each unit consisting of spins with equal magni- 
tudes. From the theoretical viewpoint, the UDC has been rigorously treated to some extent 
in the absence of distortion. ^'^^ The distorted version of this model has also been investigated 
theoretically.^^"^^^ Surprisingly, it is found that the natural mineral azurite consists of dis- 
torted diamond chains with spin- 1/2 and the magnetic properties of this material have been 
experimentally studied in detail. ^^'^^^ Other materials with the same structure have also been 
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Fig. 1. Structure of mixed diamond chain with S = 1 and r^^^ = r'^^ =1/2 



reported. ^^'^''^ 

In the present work, we investigate the thermal properties of a mixed diamond chain 
(MDC) with spins 1 and 1/2 depicted in Fig. 1. The MDC has been introduced in a previous 
paper. All the eigenstates of this model are represented in terms of the eigenstates of odd- 
length spin-1 antiferromagnetic Heisenberg chains (AFHl's) and singlet dimers in between. 
The ground state of this model has been rigorously determined in ref. 20, where it is shown 
that a series of quantum phase transitions take place among phases with different periodicities 
with spontaneous translational symmetry breakdown. In the present paper, we construct the 
full thermodynamics of the present model in terms of the eigenstates of the AFHl's with 
arbitrary odd lengths. In an appropriate parameter regime, only the eigenstates of short 
AFHl's have a dominant contribution. Therefore, using the numerical diagonalization data 
for AFHl's, we can obtain reliable numerical results for the thermodynamic properties of 
infinite mixed diamond chains. As a result, unique magnetic behaviors are predicted at low 
temperatures such as the residual entropy and Curie constant which vary from one phase to 
another. Furthermore, at phase boundaries, these quantities show different behaviors from 
those in the phases on both sides of the boundary. At higher temperatures, more exotic 
behaviors are further predicted such as the multipeak structure in the A-dependence of entropy, 
the nonmonotonic temperature dependence of the Curie constant, the multipeak structures of 
the temperature dependence of specific heat and the plateau-like behavior in the temperature 
dependence of entropy. Physical interpretations of these behaviors is also presented. 

This paper is organized as follows. In the next section, the model Hamiltonian is presented 
and the structure of the eigenstates is explained. The ground state properties are summarized 
in §3. Finite temperature statistical mechanics is formulated in §4. Low-temperature properties 
are discussed analytically in §5. The numerical analysis of physical properties such as magnetic 
susceptibility, specific heat and entropy is presented and discussed in §6. The last section is 
devoted to the summary and discussion. 
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2. Hamiltonian and Cluster Eigenstates 

The Hamiltonian of the MDC in Fig. 1 is represented as 

n = f: J(5.rf^) + S.rf ) + r«5.,, + r^S,,,) + Kr«rf \ (2.1) 

1=1 

where Si and t[°^ (a = 1, 2) are spin operators in the Zth unit ceh, J and K are the exchange 
energies, and L is the number of unit cehs. In this paper, we only consider the case in which J 
is positive, and the magnitudes of Si and r^"^^ are 1 and 1/2, respectively. The ratio A = K/J 
of the exchange energies controls the strength of frustration. Below, we set the energy unit 
J= 1. 

By defining the composite spin T/ = t[^''+tP^ for all /, Hamiltonian (2.1) is also expressed 

as 



SiTi + TiSi+i + ^Tf-l 



(2.2) 



1=1 

Then, Tf {I = 1,2,...,L) commutes with n, and T; defined by = Ti{Ti + 1) IS a good 
quantum number that takes or 1. Therefore, each eigenstate of the Hamiltonian belongs to 
one of the subspaces, each specified by {Ti}. The set {T;} is a sequence of O's and I's with the 
following structure: 

...0 1_^0 1_^0..., (2.3) 

where is the number of I's bounded by two O's. If T; = so that r|^'' and form a 
singlet dimer state, then we call the spin pair a dimer. The spins on the left of the dimer and 
those on the right of it are rigorously decoupled. Therefore, the state of the total spin chain 
with a specific set {T;} is decomposed into dimer states and local spin states between dimers 
in a tensor product form. We call such a local spin state a spin cluster state, and call the 
assembly of the spins forming the spin cluster state the spin cluster. Specifically, we call the 
spin cluster cluster-n if the spin cluster state includes n triplet-pair states (T; = 1). The series 
of cluster-n's is illustrated in Fig. 2. In particular, two successive O's imply a single spin-1 site 
bounded by them. This spin-1 site is a cluster-0, which is also called a monomer. 

A cluster-n is equivalent to an AFHl with length 2n + 1, since it consists of n spin-l's 
for T;'s and n + 1 spin-l's for Sj's, as was argued above. Thus, we can decompose the total 
spin chain into AFHl's with lengths 2?ij + 1 with ?ij > in each subspace of fixed {T;}. The 
lowest-energy state of a cluster-n-'^*) has total spin 1 from the Lieb-Mattis theorem. -^^^ Thus, 
an eigenstate of total Hamiltonian (2.2) is specified by a set of quantum numbers {rij, Oj} that 
satisfy 

Y,{ni + l) = L, (2.4) 

i=l 
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Fig. 2. Series of spin clusters bounded by two singlet dimers. A filled (open) circle represents a spin- 
1/2 (spin-1) site. An unshaded (shaded) oval represents a singlet (triplet) spin pair. A spin cluster 
is represented by a thick rectangle including triplet pairs and the spin-1 sites connected to them. 
Cluster-rt means a spin cluster including n triplet pairs. Each spin cluster has total spin 1 in its 
lowest-energy state. 



where the total number of cluster-n's is denoted by and Ui specifies the eigenstate of 
the cluster-nj. In the present work, we only consider the thermodynamic quantities in the 
thermodynamic limit. Therefore, we assume the open boundary condition for simplicity. 

The energy eigenvalue E^°^ of total Hamiltonian (2.2) is decomposed into contributions of 
clusters as 

E*°*({ni,aa,A) = J^S(n„Qi,A), (2.5) 

i=l 

where E{n, a. A) is the a-th energy eigenvalue of a cluster-n. For convenience, we have included 
the ground state energy of a neighbouring dimer — 3A/4 in E[n, a, A). This can be expressed 
using the a-th eigenvalue E{2n + 1; q) of an open AFHl with length 2n + 1 as 

E{n,a,X) = E{2n + l,a) + ^ -. (2.6) 

Correspondingly, the lowest-energy of a cluster-n can be expressed as 

n A 

Ecin, A) = EG{2n + 1) + _ - _, (2.7) 

where EQ{2n + 1) is the lowest-energy of AFHl with length 2n + 1 listed in Table I for 
< n < 4. 
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Table I. Ground state energies of odd- length AFHl's. 



2n+ 1 


1 


3 


5 


7 


9 


^G(2n + 1) 





-3 


-5.8302125227708 


-8.6345319827062 


-11.4329316403302 



Table II. Ground state phase boundaries. 



ni 712 


1 


1 2 


2 3 


3 oo 


Ac(ni,n2) 


3 


2.660425045542 


2.58274585704 


2.5773403291 



3. Ground States 

The ground states of Hamiltonian (2.1) are argued in ref. 20 for case of the arbitrary spin 
magnitude. In particular, some rigorous results and their proofs have been given. In the case 
of spins I and 1/2, all the ground states and quantum phase transitions are determined from 
numerical diagonalization data for finite AFHl's with odd numbers of spin-l's, and the full 
phase diagram is obtained. In this section, we summarize the ground state properties of this 
model for convenience. 

The ground state is a uniform array of cluster-n's with a common value of n and dimers 
in between. This state is called the dimer-cluster-n (DCn) state. The DCO state is also called 
the dimer-monomer state. The DCoo state is the Haldane state. The value of n is determined 
so as to minimize the total energy for each A. This ground state has a spatial periodicity of 
n + 1 owing to the {n + l)-fold spontaneous translational symmetry breakdown. The ground 
state energy of this DCn phase is given by 

E'^\n,X) = ^EGin,X). (3.1) 
n + 1 

Accordingly, the ground state energy in the DCn phase per unit cell is given by 

f u 1 p / u ^G(2n + 1)-A , A 

eG{n,\) = ^—EG{n,X) = — + -. 3.2) 

n + 1 n + 1 4 

If T; = 1 for all diagonal bonds, the ground state energy in the Haldane phase per unit cell is 

given by 

eG(oo,A) =2eG(oo) + ^, (3.3) 

where edoo) ~ -1.401484038971^^^ is the ground state energy of an infinite AFHl per unit 
ceU. 

By numerical analysis in our previous study,^"^ we have found that the DCO, DCl, DC2 
and DC3 phases appear successively with decreasing A from +00. The phase transition between 
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the DC(n — 1) and DCn phases takes place at A = Xc{n — 1, n) given by 

Ac(n- l,n) = {n+l)EG{2n- 1) - n^G(2n+ 1), (3.4) 

which is the solution oi eQ{n— 1, A) = eQ{n, A). The critical value is expressed by the ground 
state energies of odd-length AFHl's. On the other hand, the DCn phases with ?i > 4 do 
not appear with further decrease in A, and the DCS phase changes directly into the Haldane 
phase. The phase boundary between the DCS and Haldane phases is given by 

Ac(S, oo) = -Eg(7) - 8eG(oo) ~ 2.577S40S291, (S.5) 

which is the solution of eG(S, A) = eG(oo, A). These critical values of A are listed in Table II. 

4. Formulation of Statistical Mechanics 

4-1 Grand partition function 

We consider a state that consists of cluster-n's with lengths ni,n2,n3, • • • and njv^ from 
left to right. For convenience in calculation, we allow nj to vary freely. Instead, we introduce 
the chemical potential conjugate to the number of unit cells in order to fix the expectation 
value of the total number of unit cells as eq. (2.4). We also include the external magnetic field 
H to calculate the magnetic susceptibility x- 

The grand partition function 'E{l3, n, H) is given by 

oo A^c 

Nc=l{nj}i=l 
oo / oo \ Nc 

Nc=l \n,=0 J 
oo oo oo 

where stands for ... and Zn{P,H) is the partition function of a single 

{rij} ni=0n2=0 11^^=0 

cluster-n defined by 

Z„(/3, H)=J2 (4.2) 

a 

with the magnetization M(n, a) in the a-th eigenstate. Equation (4.1) is then represented as 

m,,H) = -^4^l^ (4.S) 
1 - =ci(/3,/U,i7) 

using the grand partition function Hci(/3, /i, ff) of a single cluster defined as 

oo 

E,i{(3,fi,H) = Y,Zn{f3,H)e^>^^^+'l (4.4) 

n=0 

Suppose that Q is an extensive physical quantity decomposed as 

Q = '^Qc\{ni,ai), (4.5) 

i=l 
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then, the grand canonical expectation value {Q) is given by 



-,00 / iVc \ 

^ Nc=l {uj} {oj} \l=l ) 
X ^~fi{E(ni,a,,\)-^{ni+\)-M{ni,Cii)H) 



00 

A^cHci(/3,/x,ii')^MQei; 



(4.6) 



iVc = l 

where (Qd) is the grand canonical expectation value of Qci in a cluster defined by 



(^c) = ^ = -r^- (4-8) 



^ 00 

""'^^ n=0 a 

Taking Nc = ^ have Qd = 1- Then, the expectation value (Nc) is given as 

oo 

c = l 

oo 

Using {Nc), eq. (4.6) is written as 

(Q) = (AT,) (Q,i) . (4.9) 

To characterize the size of a cluster-n, it is convenient to use L^i = n + 1. In fact, the sum of 
Lci's of all the clusters included in the whole chain is the total number of unit cells as shown 
in eq. (2.4). Taking L^i for Q^i in eq. (4.7), we have 

(L,i) = ^ ^(^+ l)^e-/3(i^(n,",A)-Mn+l)-Af(n,a)//)_ ^^ -^^^ 

n=0 a 

Thus, eq. (4.9) for Q = L yields 



(L) = {Nc) {Lci) = (4.11) 

1 - "Cl 

We define {{Qd)) as (Qd) / (-^d) for convenience. Equations (4.9) and (4.11) show that {{Qd)) 
is equal to the expectation value per site {Q) / {L) for the extensive quantity Q. Thus, we 
have 

Here, we determine the chemical potential fi in the thermodynamic limit. While {L) is a 
macroscopic quantity, (id) remains always finite at finite temperatures because the entropy 
gain of O(lnL) overcomes the finite energy cost in making a cut in a cluster with length L by 
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virtue of one-dimensionality. Hence, from eq. (4.11), we must require 

He, = 1 - ^ ^ 1 (4.13) 
in the thermodynamic limit. This condition fixes //. 



^.2 Expression in terms of finite-length AFHl 



We can express our formulae using the quantities associated with AFHl. For example, the 
partition function of a cluster-n is expressed as 

Z„(/3, H) = e-("-3)/3V4 z[p, H; 2n + 1), (4.14) 

where 

Z{/3, H;2n + 1) = Y^ exp[-/3(^(2n + 1, a) - HM{n, a))] (4.15) 

a 

is the partition function of the AFHl with length 2n+l in the magnetic field H. Correspond- 
ingly, the cluster grand partition function Sci(/3, /i, .ff) is expressed as 

oo 

Hci(/3, fi,H) = J2 Zn{/3, H) e^^("+i) = e^^ Hodd(/3, ^, H), (4.16) 

n=0 

where ji = n — X/4: and Hodd(/3) A) is the grand partition function of the odd-length AFHl 
given by 

oo 

Hodd(/3, f^,H) = J2 H; 2n + 1) e^f^^^+'K (4.17) 

n=0 

Therefore the grand partition function of the whole chain is expressed as 

m,,H) = /^^i^^ifML. ,4.18) 

The chemical potential jl is determined by putting = 1 in eq. (4.16) as 

X = -TlnEoddW,fi,H). (4.19) 
The cluster expectation value of the physical quantity Q can be written as 

oo 

J2 Yl «) e-^(^(2„,+l,a)-A(n+l)-Af(„„a)H) 

(Qci) = . (4.20) 

^ ^ g-/3(S(2n+l,a)-A(n+l)-Af(n,a)H) 
n=0 a 

This expression is useful in numerical calculation, because once the full spectrum of the finite 
length AFHl is calculated numerically, we can evaluate the physical quantities of model (2.1) 
using eqs. (4.12), (4.19) and (4.20) without further massive numerical calculations. 
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4.3 Magnetic susceptibility 

The expectation value of the magnetization M in the magnetic field H is obtained by 
setting Q = M in eq. (4.12) as 

(M) = (L) ((M,i)) = (L) (4.21) 

where {Md) is given by eq. (4.7), using M^i for Q^i- Within the first order in H, the in- 
dependence of /U can be neglected, because fi is an even function of H from symmetry consid- 
eration. Therefore, we find 

00 

lim = hm 



-l3{E{n,a,\)-fi{n+l)) 

n=0 a 

(Ml) 



(4.22) 



where the expectation value {M"^^ is taken in the absence of the magnetic field. Consequently, 
the magnetic susceptibility x of total chain is given by 

X= HmM = i^((A42,)). (4.23) 

4-4 Entropy 

Here and hereafter, we take H = 0. The entropy S is evaluated using the formula 



S = T^TlnH 
dT 



= ln(iVe) + ^[(i^ci)-^(Lci)]. (4.24) 



In the thermodynamic limit, In (Nc) can be neglected in comparison with (A'c). Therefore, we 
find 

5 = ^(((i?ci))-M)- (4.25) 

4.5 Internal energy and specific Heat 

The expression for the internal energy U of the total chain is obtained by setting Q{n, a) = 
E{n,a,X) in eq. (4.12) as 

n a 

Differentiating eq. (4.26) with respect to T with fixed L, we have the following expression 
for specific heat per unit cell: 



C _ J_ 

T ~r2 



{{El)) -{{Ed)) {{EM) 
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dp 



Since the differentiation of eq. (4.13) for H = with respect to /3 yields 



9/3 



{{Ed)), 



(4.27) 



(4.28) 



1 



[{{E^,)) - 2{{E,,)){{EM) + {{Eci)){{Ll))] . 



(4.29) 



eq. (4.27) is written as 

C _ 

5. Low- Temperature Limit 

5.1 Low-temperature behavior of each phase 

The lowest-energy state of each cluster-n has total spin 1 owing to the Lieb-Mattis the- 
orem. Hence, L/{n + 1) spin-1 degrees of freedom remain in the DCn phase at T = 0. 
Accordingly, we respectively have the magnetic susceptibility and residual entropy in the 
DCn phase as 

2 2 



X^^{Nc 



5~ (ln3 



3{n+l)T 
ln3 



{L) 



n + l 



(L). 



(5.1) 



(5.2) 



Equation (5.1) means that the Curie constant in each DCn phase is fixed to an n-dependent 
value (2/3)L/(n -|- 1). Therefore, we have a stepwise A-dependence of xT, as shown by the 
solid line in Fig. 3. 

5.2 Low-temperature limit on the phase boundary 
5.2.1 Magnetic susceptibility 

Near the phase boundary between the DC(n — 1) and DCn phases, only the lowest-energy 
states of cluster-(n — 1) and cluster-n contribute to Hodd- Therefore, we find 



Sodd(/3, /i) = [^(/3; 2n - 1) + Z(/3; 2n + 1)6^^^^] e^'^^ 

where 

Z(^;2n± 1) = 3exp(-/3£'G(2n± 1)). 
According to the formula (4.19), the chemical potential /i satisfies the relation 



3(1 + ^-P(EG(2n+l)-EG(2n-l)-fL) 



A = EG{2n - 1) - /in - Tin 



On the phase boundary, we combine eqs. (3.4) and (5.5) to find 



X 

3 



1 — X 



where 



1 + exp(-/3(£;G(2n + 1) - EG{2n - 1) - /i) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 
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Table III. x, x^/L and S / L ai \ = \c{n ~ l,n) in the low-temperature limit. 



n — 1 n 


X 


XT/L 


S/L 


1 

1 2 

2 3 


0.791287847477920 
0.677814645373914 
0.627507207887062 


0.55155122356933 
0.28708589748815 
0.19767771430847 


1.33270576282026 
0.74374685074522 
0.52153747001256 



which implies that 



ii = T\J—^+EG{2n + l)-EG{2n-l). (5.^ 



The solutions of eq. (5.6) are tabulated in the second column of Table III. For A = Ac(0, 1), 
this equation gives an analytic value of a; = (\/2T — 3)/2 (= 0.791). 

Using the values of x in Table III, the low-temperature susceptibility on the phase bound- 
ary between the DC(n — 1) and DC?i phases is expressed as 

(L) o 2(L) e-^EG(2n-l) , -/3(£G(2n+l)-A) 



T " 3r ^e-/3£;G(2n-l) _^ l)g-/3(£;G{2n+l)-/i) 

For A = Ac(0, 1), this equation gives an analytic value of xT jh = {\/21 + 7)/21 (= 0.552). in 
the low-temperature limit. The numerical values of x^/L for all A = Ac(n — 1, n) are tabulated 
in the third column of Table III. 

5. 2. 2 Entropy 

The entropy S is also expressed by x given by eq. (5.6). Using eqs. (4.25) and (5.8), S is 
calculated as 



E^in - 1, A)e-^^G(2n-l) ^ ^^/^^ ^y-p^EG{2n+l)--^) 

^ 



ne 



-PEG{2n-l) + (n + l)e-/3(-BG(2n+l)-A) 



_ (L) (n + l)EG{2n - 1) - nEoi'^n + 1) - A - ^ 

T n+l-x ^ ' X ' ^ ' ' 

Since the first term vanishes as indicated in eq. (3.4), we finally obtain 

S = -{L)ln^—^. (5.11) 

X 

The numerically estimated values of 5 are tabulated in the fourth column of Table III. For 
A = Ac(0, 1), this equation gives an analytic value of S/L = \n{{\/21 + 3)/2} (= 1.333). Note 
that the entropy at A = Ac(n — l,n) is larger than those in the DCn and DC(n — 1) phases 
because of the excess entropy due to the mixing of two types of clusters. These results can be 
reproduced by the direct combinatory calculation in the ground state, as shown in Appendix. 
At A = Ac (3, oo), no finite cluster-n's with n > 4 can coexist with cluster-3's. Although the 
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excitation energy of macroscopic cluster-n's with n ~ 0{L) vanishes in the thermodynamic 
Hmit, the introduction of a macroscopic cluster reduces the macroscopic amount of entropy. 
Therefore, no macroscopic clusters can appear at A = Ac(3,oo) and residual entropy remains 
the same as that in the DC3 phase. 

6. Numerical Calculation 

6. 1 Estimation of contribution from large clusters 

We numerically calculate physical quantities using the general formulae derived above. 
Owing to the computer memory limitation, we can only include the contribution of cluster- 
n's where the cluster size n is not very large. In the following calculations, we denote the 
largest cluster size Ti employed in the numerical calculation as nj^^x;* 

We examine the missing contribution of discarded clusters larger than rimax by estimating 
the entropy in the high-temperature limit. From eq. (4.25), the entropy per unit cell is 
expressed as S/L = /3{{Eci)) — Pfi. The first term (3{{Eci)) vanishes in the limit of T — )• oo, since 
the physical quantity {{Ed)) is finite. Thus, denoting the limiting value of S/L as Soo/L, we 
have 

— = - lim (3n. (6.1) 

L /3^0 

We use this relation to estimate the approximate entropy in the high-temperature limit formed 
by cluster-n's with 1 < n < nmax- We denote it as 5oo(nmax)i and then 5oo(oo) is the exact 
value of entropy in the high-temperature limit. In this limit, we have Zn = 1 = 3^"^^ 
from eq. (4.2). Hence, for nmax; the cluster grand partition function is 

"max 1 _ /'Q„/9/l\"max + l 

= y 32"+ie^^("+i) = 3e^^ ^ \ . (6.2) 

n=0 

From eq. (4.13), we require Hd = 1. Then eq. (6.2) reduces to the following equation: 



1 + I {9e^^ 



(6.3) 



12 

The solution of this equation provides the estimation iS'oo(nmax)/-^^ for each nmax- 

The 

analytic solutions for n^ax = 0, 1 and 2 are 5'oo(nmax)/-^ = ln3, ln[(3/2)(\/T3 + 1)] and 21n3, 
respectively. Solutions for nmax > 3 are numerically obtained. Also, for nmax oo, eq. (6.3) 
gives Sooioo)/L = ln3 + ln4, which agrees with the entropy calculated from the number of 
states 12 for one spin-1 and two spin-l/2's per unit cell. 

We present results for 5'oo(nniax)/-^ obtained from eqs. (6.1) and (6.3) in Table IV. The 
second column contains the values of >S'oo(nniax)/-^ for 1 < n^ax ^ 5. The third column 
contains those of the deviations A5oo(nmax)/-^^ = 5'oo(oo)/L — Scx]{nraa.x)/ L. For comparison, 
we show the exact value of the high-temperature entropy in the last row. This table shows 
that the missing entropy is only about 3% for nmax = 5, even in the high-temperature limit. 
Therefore, we satisfy ourselves with nmax = 5 and use only results for nmax = 4 for comparison. 
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Table IV. Contribution of only cluster-n's with < n < rimax to the entropy in the high-temperature 
limit. The second column contains the approximate entropy per unit cell for each rimax- The third 
column contains the deviation from Sao{oo)/L = ln3 + ln4, which is also shown in the last row. 



'^max 


5'oo(^max)/ L 


A5oo(nmax)/-^^ 





1.098612 


1.386294 


1 


1.932727 


0.552179 


2 


2.197225 


0.287682 


3 


2.315806 


0.169100 


4 


2.378360 


0.106547 


5 


2.414692 


0.070215 


oo 


2.484907 






if necessary. 

6.2 \- dependence of magnetic susceptibility and entropy 

Figure 3 shows the A-dependence of the magnetic susceptibility times the temperature 
at various temperatures estimated with rimax = 4 and 5. This quantity approaches the Curie 
constant in the low-temperature limit. The stepwise dependence of xT in the low-temperature 
limit is smeared out by thermal fluctuation at finite temperatures. For A < Ac(3,oo), the 
ground state is the Haldane phase with n — )■ oo. Therefore, the present approximation is not 
reliable. Actually, the result strongly depends on nmax(4 or 5) in this region. In contrast, our 
approximation is fairly reliable for A > Ac(3, oo), where the ground state is an assembly of 
finite-size clusters. In this region, the excited states with large values of n are unfavorable not 
only energetically but also entropically, because the number of such excited states is small, as 
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Fig. 4. A-dependence of entropy with Timax = 5 for Ac(3, oo) < A < 3.1 at various temperatures. 

estimated in §6.1. We actually find that the results in this region are insensitive to "/T/max; 
shown in Fig. 3. Therefore, we concentrate on the region A > Ac (3, oo) and fix /^max 

at 5. 

The A-dependence of entropy is shown in Fig. 4. The peaks due to the excess residual 
entropy at the phase boundaries are smeared out by the thermal fluctuation with an increase 
in temperature. 

6.3 Cluster excitation energy 

At finite temperatures, spin clusters other than the cluster-n are excited in the DCn 
phase. One of the simplest excitations is a cluster-n' with n' ^ n in its lowest-energy state. If 




Fig. 5. A-dcpcndcncc of the excitation energy AEG{n' ;n) of cluster-n' in the DCn phase, where 
6X = X — Ac(3,cx)). The size n' of the excited cluster is indicated for each curve for < n' < 5. 
The line AEq = J (=1) is plotted by a dotted line as a guide for the eye. The triangles on the 
abscissa indicate the values of the ground state phase boundaries, and each ground state phase is 
indicated under the abscissa. 
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a cluster-n' is created in the DCn ground state, the system size L changes by n' — n. Hence, the 
excitation energy with fixed L is given by subtracting the contribution of chemical potential 
as lS.EQ{n'; n) = EQ(n', A) — EQ(n, A) — (n' — In the low-temperature limit, the condition 
(4.13), which determines fi, reduces to the form 

-^oo,n)= lim 3e-/3(-^5G(n,A)-M(n+i)) = (5.4) 

which implies that /x = EQ{n)/{n + 1). Thus, the excitation energy required to create a 
cluster-n' in its lowest-energy state is written as 

r ' ^ , , , / ^G(2n' + 1)-A ^G(2n + 1)-A | 

AEcin ; n) = n + 1) <^ — — S . 6.5) 

I n + I n + 1 I 

The A-dependence of A.EQ,{n' ]n) for various values of n' is shown in Fig. 5, where we use 
5X = \ — Ac(3, 00) instead of A itself. Although this quantity is simply linear in A, we employ 
the log- log plot to magnify small energy differences. 

6.4 Temperature dependence of magnetic susceptibility 

The temperature dependence of xT is presented in Fig. 6 for various values of A. It shows a 
distinct nonmonotonic behavior. By calculating xT with rimax = 3, 4 and 5, we have confirmed 
that the behavior is insensitive to the choice of fimax- Therefore, this is not an artifact of the 
present approximation. This behavior is understood by comparing Fig. 6 with Fig. 5 as follows. 

We investigate xT near each phase boundary Ac(n — l,n). If A approaches Ac(ri — l,n) 
in the DCn phase (A < Xc{n — l,n)), the excitation energy AEq^h — l;n) of cluster-(?i — 1) 
comes down to zero, as shown in Fig. 5. Hence, cluster-(n — l)'s mix with cluster-n's at low 
temperatures. A cluster-(n — 1) has the same spin unity as a cluster-n, and the cluster length 




Fig. 6. Temperature dependence of xT / L within the phases (open squares) and on the phase bound- 
aries (filled squares) with nmax = 5. 
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of a cluster- (n — 1) is shorter by a unit. Accordingly, such excitations increase the density 
of alive spins, which lead to an increase in at low temperatures. This behavior is seen, 
e.g., for A = 2.9 near Ac(0, 1) = 3 in 10"^ < T < 10^^ in Fig. 6. On the other hand, if A 
approaches Ac(n — in the DC(n — 1) phase (A > Ac(n — 1,?^)), the excitation energy 
A£'g(?t-; ^^— 1) of cluster-n comes down to zero, as shown in Fig. 5. Hence cluster-n's mix with 
cluster-(n— l)'s at low temperatures. Such excitations decrease the density of alive spins, and 
lead to the decrease in xT- This behavior is seen, e.g., for A = 2.7 near Ac(l,2) = 2.660 in 
10"^ < T < 0.5 X 10^1 in Fig. 6. 

On the phase boundary of A = Ac(n — the energies per unit cell of cluster-(n — 1) 

and cluster-n degenerate. In this case, the next-lowest excitation controls the temperature 
dependence of xT- At A = Xc{n,n — 1), the candidates for the lowest excitation energy are 
AEq^ti — 2;n) and AEQ{n + Using eqs. (3.4) and (6.5), these are written as 

AEcin - 2; n) = ^G(2n + 1) + £^G(2n - 3) - 2^G(2n - 1), (6.6) 

AEcin + 1; ii) = ^G(2n + 3) + ^G(2n - 1) - 2^G(2n + 1). (6.7) 

Because of their second difference forms, they are independent of the bulk term proportional to 
n and the constant boundary term in EQ{2n+l). Accordingly, A£^g('^~2; n) and Aii^G(^+li n) 
are only determined by the n-dependent boundary term resulting from the interaction between 
the boundary spins of the cluster-n. Since this interaction decreases with n, A£'g(^ + li ?^) is 
lower than A£'g(^ — 2; n). Figure 5 actually shows that, at each phase boundary Ac(n — 1, n), 
the next-lowest cluster is cluster- (n -|- 1) with an excitation energy of Aii^G(^ + Ij^T-)- The 
mixing of cluster-(n + l)'s into cluster-(n — l)'s and cluster-ji's decreases the density of alive 
spins. Therefore, even at the phase boundary, xT decreases with T at low temperatures. 

At higher temperatures, clusters with higher energies affect x^- Figure 5 shows that 
cluster-(n — 2)'s come into play for A > Ac(n — l,n), e.g. an excitation of n' = appears in 
DCl (A > Ac(l, 2)) phase and at the phase boundary A = Ac(l, 2). The mixing of these shorter 
clusters increases x^, which is decreased by cluster-(n + l)'s at low temperatures. 

At a general temperature less than J, short and long clusters with various excitation ener- 
gies contribute to x^, and the shorter ones increase and the longer ones decrease xT- Which 
clusters mainly contribute to x"^ depends on temperature, and this qualitatively explains the 
nonmonotonic temperature dependence of xT. For T > J, xT increases monotonically toward 
the free-spin limit x^/L = 2/3 + 1/2 = 7/6. The temperature dependence of xT is insensitive 
to A in this regime. 

6. 5 Temperature dependence of specific heat and entropy 

The temperature dependence of specific heat is presented in Fig. 7. In general, specific heat 
shows a multipeak structure. This structure is also understood on the basis of the excitation 
spectrum AE'g (n';n). 
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Fig. 7. Temperature dependence of specific heat (a) in each phase and (b) on each phase boundary 
by calculation with rimax = 5. 




Fig. 8. A-dependence of the peak temperature Tpcak of specific heat where SX = X — Xc{3,oo). The 
size n' of the cluster responsible for each peak is indicated for each curve with rimax = 5. The 
ground state phases are indicated below the abscissa. The triangles are the ground state phase 
boundaries. 
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Fig. 9. Temperature dependence of specific heat around the phase boundary of Ac (0,1) = 3 with 
"max = 5 for (a) A > 3 and (b) A < 3. 



The A-dependence of the peak position is shown in Fig. 8. The uppermost peak is common 
to ah values of A, and originates from the bulk excitations of cluster-n's, which are equivalent 
to finite AFHl's. The positions of the lower peaks strongly depend on A. When A approaches 
a phase boundary from any side, a peak approaches T = 0, as shown in Fig. 8. By comparing 
Fig. 8 with Fig. 5, such a peak is interpreted to originate from the excitation of a cluster- 
(n + 1) or cluster-(n — 1) with an excitation energy of AEq^u it 1; n) in the DCn ground state. 
Approaching the phase boundary, AEQ{n + l;n) comes down to zero for A > Xc{n,n + 1) 
and AEQ{n — l;n) comes down to zero for A < Xc{n — l,n). Hence, at the phase boundary, 
cluster-n's and cluster-(n — l)'s mix, and yield the excess residual entropy at T = 0. We show 
the temperature dependence of specific heat around the phase boundary of Ac(0, 1) = 3 in 
Fig. 9, where specific heats are plotted at A = 2.99 to 3.008 with an interval of 0.002. The 
lower peak shifts toward T = as A approaches Ac(0, 1) = 3 from both above and below. At 
A = 3, the lower peak vanishes. Also, in the neighbourhood of other phase boundaries, specific 
heat behaves similarly. 

Entropy is estimated using the formula (4.25). The temperature dependence of the entropy 
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Fig. 10. Temperature dependence of entropy around (a) A — Ac(0,l), (b) A = Ac(l,2) and (c) 
A = Ac (2, 3) with rimax = 5. The open circles, filled squares and open squares represent the 
values for A = Ac, A > Ac and A < Ac, respectively. The dotted line indicates the correct value of 
entropy in the high-temperature limit. The solid lines are the values of residual entropy in each 
phase. The broken lines are the values of residual entropy on the phase boundaries. 

is shown in Fig. 10, where panels (a), (b) and (c) are for regimes around A = Ac(0, 1), A = 
Ac(l, 2) and A = Ac (2, 3), respectively. We find that entropy has a plateau- like behavior around 
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the value of the residual entropy on phase boundaries. This behavior is explained as follows. 

Between the temperatures corresponding to J and AEq{0; 1) (A < Ac(0, 1)) or AEq{1;0) 
(A > Ac(0, 1)), there exists a temperature regime where bulk excitations of the finite-length 
AFHl are not excited and several cluster-n's with n > 0, which have lower excitation energies 
than the Haldane gap, are still equally excited in their lowest-energy states. In this regime, the 
entropy is close to the residual entropy at A = Ac(0, 1) because the contribution of cluster- Ji's 
with n > 2 is small. For large A ^ Ac(l,2), this behavior continues until the true ground 
state is chosen at lower temperatures, as shown in Fig. 10(a). With a decrease in A, however, 
AEq(0; 1) increases while AEq^h; 1)'s with n > 2 decrease. Therefore, there exists a tem- 
perature regime where cluster-n's with n > 1, which have lower AEq^h; 1) than A£'g(0; 1), 
are still equally excited in their lowest-energy states. The entropy of this state is close to the 
residual entropy at A = Ac(l,2), as shown in Fig. 10(b), because the contribution of cluster- 
n's with n > 3 is small. Similarly, the temperature dependence of entropy has plateau-like 
structures around the value of the residual entropy on the phase boundary at Ac (A; — 1, A;), if 
A ~ Ac(n — 1, n) with n > k. 

7. Summary and Discussion 

The statistical mechanics of the MDC with spins 1 and 1/2 is formulated rigorously. 
The low-temperature behaviors of the magnetic susceptibility, specific heat and entropy are 
analytically investigated. The Curie constant and residual entropy vary depending on the 
ground state phase. At finite temperatures, these quantities are obtained using the exact 
numerical diagonalization data for odd-length AFHl's. The finite-temperature behavior of 
physical quantities also varies reflecting the ground state properties and excitation spectra. 
Magnetic susceptibility shows a nonmonotonic temperature dependence. Specific heat has a 
multipeak structure. At low temperatures, entropy as a function of A has peaks around the 
ground state phase boundaries refiecting the ground state excess entropy. Entropy as a function 
of temperature exhibits plateau- like structures depending on A. An intuitive explanation for 
these exotic features based on the cluster picture is given. Similar behaviors are also found in 
the frustrated Ising model on a diamond hierarchical lattice. 

Although no real material described by the MDC model has been known so far, the MDC 
with distortion will hopefully be realized. It is well known that the natural mineral azurite, 
in which Cu ions carry spin- 1/2 degrees of freedom, is a distorted version of the spin- 1/2 
UDC. Another material for the distorted UDC with spin-1/2 is also known. These facts 
encourage the search for materials described by the distorted version of the MDC. As long 
as the distortion is weak, the effect of distortion is smeared out at finite temperatures and 
the phenomena predicted in the present work would be observable. On the other hand, the 
ground state properties are sensitive to lattice distortion because the DCn ground states of 
the present model have a macroscopic number of effectively free spin-1 degrees of freedom. 
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A variety of quantum phases can emerge in the presence of perturbation. Therefore, it is 
also important to investigate the effect of lattice distortion on the ground state of MDC's 
theoretically. Our preliminary study suggests the presence of a variety of exotic phases and 
phenomena, such as quantized and partial ferrimagnetic phases, Haldane phases with broken 
translational symmetry and an infinite series of quantum phase transitions depending on the 
type of distortion. Further investigation of these phenomena is in progress and will be reported 
in future publications. 

The interchain interaction is another important perturbation that controls the low- 
temperature behavior of the MDC, if it is realized as a real material. It is naturally expected 
that the spin-1 degrees of freedom, which remain paramagnetic in the ground state of a sin- 
gle chain, become ordered in the presence of an interchain interaction. The competition of 
magnetic ordering against the various phenomena reported in this work is also an interesting 
issue. 

The numerical diagonalization program is based on the package TITPACK ver.2 coded by 
H. Nishimori. The numerical computation in this work has been carried out using the facilities 
of the Supercomputer Center, Institute for Solid State Physics, University of Tokyo and the 
Supercomputing Division, Information Technology Center, University of Tokyo. This work 
is supported by a Grant-in- Aid for Scientific Research on Priority Areas, "Novel States of 
Matter Induced by Frustration" , from the Ministry of Education, Science, Sports and Culture 
of Japan, and Fund for Project Research in Toyota Technological Institute. 

Appendix: Combinatorial derivation of the residual entropy 

We reproduce the entropy in §5.2.2 by direct combinatory calculation. On the phase 
boundary between the DC(n — 1) and DCn phases, two kinds of clusters, cluster-(n — 1) and 
cluster-n, coexist. We denote the number of cluster-(n — l)'s by Nn-i and that of cluster-n's 
by Nn- Under the constraints (n — l)A^„_i + nNn = L and Nn~i + Nn = Nc, the total number 
of the allowed states is given by 



in the thermodynamic limit. Here, we have taken into account the 3-fold degeneracy of the 
spin-1 lowest-energy state for each cluster. Setting Nc = Lc, Nn^i = NcX and Nn = Nc{l — x), 



W = 




n 



'Nr,-1+N, 



n 



(A-l) 



we have c 



1 



. Then, eq. (A-l) is written as 



n + 1 — X 



W = 



•Lc 



(A-2) 
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and entropy is calculated as 

S = \nW = [xlnx + (1 - x)ln(l - x) - ln3] . (A-3) 

n + 1 — X 



Optimizing S with respect to x, we have 
dS L 



dx {n + 1 — x] 



[xlnx + (1 — x)ln(l — x) — ln3] 



^ -In^— = 0. (A-4) 



n + 1 — X 1 — X 

This equation is equivalent to eq. (5.6) determining x, and eq. (A-3) reduces to the expression 
(5.11) for residual entropy. 
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